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Outflows and mass accretion in collapsing dense cores with 
misaligned rotation axis and magnetic field 
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ABSTRACT 

Outflows and jets are intimately related to the formation of stars, and play an impor- 
tant role in redistributing mass, energy and angular momentum within the dense core 
and parent cloud. The interplay between magnetic field and rotation is responsible for 
launching these outflows, whose formation has been generally carried out for idealized 
systems where the angle a between the rotation axis and large-scale magnetic field is 
zero. Here we explore, through three-dimensional ideal magneto-hydrodynamic simu- 
lations, the effects of a non-zero a on the formation of outflows during the collapse of 
dense pre-stellar cores. We find that mass ejection is less efficient for increasing angle 
a, and that outflows are essentially suppressed for a ~ 90°. An important consequence 
is a corresponding increase of the mass accreted onto the adiabatic (first) core. In ad- 
dition, mean flow velocities tend to increase with a, and misaligned configurations 
produce clumpy, heterogeneous outflows that undergo precession, and are more prone 
to instabilities. 
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1 INTRODUCTION 

The ejection of plasma in the form of more or less collimated 
bi-polar flows (jets) is arguably one of the most spectacu- 
lar displays staged by a forming star. Over a million years, 
these outflows trace the creation of stars from their embry- 
onic emergence as proto-stars, still enveloped in their dense 
parental cloud, to their birth as T Tauri stars surrounded 
by proto-planetary discs. The existence of a close rela- 
tion between the ejection and accretion of matter through- 
out the (low-mass M* < 2Mpj) star's format i on history 
JCabrit fe BertoutJ Il992l : iHartigan et~ai1 Il995l : IWu et all 



2004) suggests the presence of a single launching mecha- 
nism, which is now widely thought to rest on the presence of 
a large-scale magnetic field mediating the extraction of grav- 
itational energy from the accreting plasma, and redirecting 
mass and energy, in the form of kinetic and Poynting flux, 
into bi-polar jets. The fundamental magneto-hydrodynamic 
(MHD) mechanism for launching and collimating an axisym- 
me tric, self-similar out f low w as detailed in analytical work 
by iBlandford fc Pavnei Jl982l^l and applied to young stellar 
systems bv lPudritz fc Normanl l|l98&t ); i t was later extended 
to include non-self similar solutions iPelletier fc Pudrit j 
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(1992) and global disc-jet solutions bv lFerreiral (|l997l ). The 
basic ideas underlying MHD jet launching were confirmed 
in time-depende nt two-dimensional, axis ymmetric numeri- 
cal simulation s (lllchida fc Shi bata 19 85[) which ei ther in- 
cluded a disc (|Kudoh et al. 19981 : Zanni et al. 2007h or pre- 
scribed it as boundary conditions (lOuyed fc Pudrit j 1 19971 ; 
Ustvugova et all 119991 ; I Anderson et all 120051 : iFendtl I2OO6I ; 
Pudritz et a.l.ll2006h . Furthermore, the early eje ction of bipo- 
lar o utflows was also demonstrated in two- (|Mellon fc Lil 
2008) and three-dimensional MHD simulations of the col- 
lapse of pre-stellar dense cores (|Hennebelle fc Fromand 
2008), which notably included the se cond collapse and the 
formation of jets from th e protostar |Machida et al.| [2006b ; 
iBaneriee fc Pudritdl200d ; iMachida et al.ll2008r i. In fact, far 
from being a passive tracers of star formation, these out- 
flows play an active role which is still to be fully understood. 
Examples include t he removal of angula r momentum from 
the accreting flow jBacciotti et al.l 120021 '). determining the 
efficien cy of star formation by d ispersing the accreting en- 
velope (jMatzner fc McKe c 2000), and sustaining turbulence 
in their local environment ( Nakamura fc Lill2007l ). 



In general, the modelling of collapsing cores and jet 
launching has been carried out under the simplifying as- 
sumption that the angle a between the rotation axis and 
the large-scale magnetic field, is zero. Exceptions are the 
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simulations of iMachida et alj i|2006al ); |p"rice fc Bate! (|2007h : 
iHennebelle fc Ciardil l|2009l ), which however focused mostly 
on the core and disc dynamics. Here we present simulations 
of the collapse of pre-stellar dense cores with misaligned ini- 
tial configurations (a — 0° — 90°), and show that the angle 
q is fundamental in determining not only the properties of 
the outflow, but also the mass accretion onto the core. In 
particular we observe the gradual decrease in the mass ejec- 
tion with increasing angle a, until its total suppression for 
nearly perpendicular configurations. 



2 RESULTS 

2.1 The numerical model 

We follow numerically the gravitational collapse of dense 
pre-stellar cores, up to, and including the formation of the 
first (or adiabatic) cor e. The numeric a l simulations were per - 
formed with Ramses |Tevssieij |2002| ; iFromang et all [2006). 
an adaptive mesh refinement code which uses a Godunov- 
type scheme and constrained transport to solve the ideal 
MHD equations. Throughout the simulations the Jeans 
length is resolved with at least 10 cells, and a HLLD solver is 
employed. The initial conditions consist of a one solar mass 
core whose density profile resembles the observed cores, and 
is given by n(r) = n c / [l + ( r / r o) 2 ] where ro ~ 1000 AU is 
the inner cloud radius. This isothermal cloud (T ~ 10 K) is 
placed inside a warm and diffuse medium in pressure equi- 
librium with the cloud edge, with a contrast of 10 between 
central, n c = 8 x 10 6 cm -3 , and edge densities. The cloud 
is initially in solid body rotation and threaded by a uniform 
magnetic field along the z-axis, whose intensity is propor- 
tional to the column density of the cloud. The rotation axis 
is in the x — z plane and makes an angle a with respect 
to the magnetic field. For ease of presentation we define, 
in addition to the Cartesian coordinates (x, y, z), cylindrical 
coordinates (zn, <j>, Z) with Z parallel to the initial rotation 
axis. All simulations are characterized by four parameters: 
the ratio of rotational over gravitational energy (~ 0.03), the 
ratio of thermal over gravitational energy (~ 0.25), the de- 
gree of magnetization /i, and the angle a. Only the effects of 
different /x, the mass-to-flux over critical mass-to-flux ratio, 
and a are investigated in this work. 



2.2 Outflow dynamics, precession and instabilities 

The general outflow formation and dynamics are pre- 
sented in Figure Q] As the initially spherical, magnetized 
pre-stellar core undergoes gravitational collapse it flattens 
preferentially along the magnetic field lines, producing a 
dynamically-collapsing, magnetized disc-like structure or 
pseudo-disc. Following the isothermal collapse phase, an adi- 
abatic core with densities > 10 10 cm -3 and a radius ~ 10—20 
AU, forms at the centre of the infalling envelope. The en- 
suing build-up of a centrifugally supported disc (or simply 
disc), with characteristic diameters ~ 50 — 200 AU, depends 
on the efficiency of magnetic torque s to remove angular 
mome ntum from the pseudo-disc. In IHennebelle fc Ciardil 
(2009) it was argued that the magnetic braking efficiency is 
proportional to hT 1 ^ 2 , where the characteristic scale-height 
h of the disc is in turn determined by the angle 01 between 



the initial magnetic field and rotation axis. Increasing their 
misalignment produces thicker pseudo-discs which are less 
efficiently braked by the magnetic field, and thus lead more 
easily to the formation of centrifugally supported discs. In 
general we find that outflows are launched for all angles 
a < 80° independently of the existence of a centrifugally 
supported disc. However when a ~ 90° the ejection of mat- 
ter is essentially suppressed, even when a disc is present. 
The remaining outward motions observed in the perpendic- 
ular configuration are particularly interesting because mass 
ejection does not produce either a magnetic cavity, or a jet. 
In fact, mass is not launched from the core, but the out- 
flow traces dense regions of the disc which are at a large 
radii from the core (~ 200 AU) and are being radially ex- 
pelled by the highly twisted magnetic field. The outflowing 
gas propagates perpendicularly to the rotation axis (see Fig- 
ure [2]) and in the plane of the disc, with radial ejection and 
azimuthal speeds, « OT « v<j,. 

For ct < 80° the formation of bipolar outflow begins 
as the magnetic field lines, connecting the outer regions 
of the cloud to its rapidly rotating inner regions, undergo 
strong shear which results in a significant azimuthal com- 
ponent (B^) of the magnetic field being generated close 
to the adiabatic core and disc (if present). The growing 
magnetic pressure gradient accelerates the plasma, inflat- 
ing bi-polar magnetic cavities within the infalling envelope. 
At this stage the ratio of MHD Poynting to kinetic flux 
magnitudes, F = v±B 2 /2-Kpv 3 , where v± is the component 
of the velocity perpendicular to the magnetic field, and the 
plasma-/?, defined as the ratio of the thermal to magnetic 
pressure, are typically r ~ 2 - 10, and /? ~ 10~ 3 - 0.1 
in the magnetic cavity. Such "magnetic tower" structure 
is similar to that described by iLvnden-Belll J2003h and re- 
produ ced in scaled laboratory experiments bv lLebedev et all 
(2005), and it is inflated by the magnetic field and is not 
mechanically-driven by wide-angle winds. At the base of the 
cavity the continuous generation of provides the Poynt- 
ing flux powering an outflow which originates from the core 
and a region extending several AU around it. Depending 
on the ejection efficiency, which is measured by the ejec- 
tion index £ = din M a /dlnzu with M a the mass accre- 
tion rate through the disc, this outflow can be described 
as either magneto-centrifu gally or magnetic pressure driven 
jFerreira fc Pelletiedll995n . As we shall see in the next sec- 
tion, the ejection efficiency increase in time and both launch- 
ing regimes are attained during the simulations. The colli- 
mation of this "disc" wind i nto a jet depends o n the radial 
distribution of currents fe.g. lPudritz et aLll2006l ). and there- 
fore on the width of the magnetic cavity itself, which takes 
up the "return current" and the magnetic stresses ass ociated 
with the Bj, component of the magnetic field (e.g. ISpruitl 
2009). However, the collimation of the magnetic tower de- 
pends in turn on the distribution of gas and magnetic field 
in the surrounding environment, which is swept up into a 
shock layer delineating the walls of the expanding magnetic 
tower. Therefore the medium through which the magnetic 
tower expands plays as well an important role on the over- 
all collimation of the outflow. Figure [2] shows the outflows' 
three-dimensional structure for different a, indicating that 
increasing the misalignment leads to the narrowing and bet- 
ter collimation of the whole magnetic tower. This is due to 
the increasing scale-height of the pseudo-disc with a, which 
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Figure 1. Slices in x — z plane of the particle number density (cm 3 on a logarithmic scale) for the case fi = 5. The first row shows for 
a = 20° the times: (a) 20300 years, (b) 21300 years, and (c) 22200 years. The second row shows at ~ 27000 years the angles: (d) a = 0°, 
(e) a = 45° and (e) a = 70° . Note the use of the different scales between images. 



provides confinement over a longer lenght of the magnetic 
tower. In addition, because the outflow is inclined with re- 
spect to the pseudo-disc, it tends to propagates into it, lead- 
ing again to a better lateral confinement, albeit asymmetric. 
The difference between the aligned and misaligned cases is 
also evident in Figure [TJi and e, where the base of the mag- 
netic cavity, which is still embedded in the pseudo-disc, is 
indeed narrower for the misaligned case. 

The angle a has also a strong impact on the homo- 
geneity of the outflows. We find that for misaligned con- 
figurations the outflows undergo precession around the z- 
axis, and that this effect is smaller for larger values of 
a, settling within ~ 2000 years to an approximately con- 
stant angle 8, which we remind is initially zero. Nutation 
is very small, and the initial angle between the magnetic 
field and rotation axis, a, remains approximately constant. 
Typical values of (a, 0), for fx = 5, are (20°, 60°), (45°, 25°), 
(70°, 20°) and (80°, 10°). In addition, while the magnetic 
tower is initially launched perpendicular to the pseudo-disc, 
thus approximately parallel to the magnetic field, the out- 
flow is ejected parallel to the rotation axis of the precessing 
core/disc, resulting in asymmetric cavities as shown in Fig- 
ure [T] and [2] At later times the spatial variability observed 
in the outflows is mainly the consequence of growing non- 
axisymmetric current-driven modes and episodic ejections. 
The unsteady behaviour appears to be related to the in- 
creasing mass ejection efficiency , and build up of B^, in the 
core/disc (I Anderson et al.ll2005l ). The precession present in 
the misaligned models further helps to seed large kink-like 
perturbations in the jet body, thus promoting the growth 
of the instability, which is otherwise artificially suppressed 
in the aligned case. The presence of a kink-like instability 



is indeed expected in jets confined by a helical magnetic 
field, and in detailed three-dimensional simulations it was 
also seen to lead to co nsiderable distort ions of the jet body 
and its fragmentation (|Moll et al.l l2008). In our simulations 
we find that the outflows are above the Kruskal-Shafranov 
stability threshold, with typical values of the magnetic pitch 
evaluated along single field lines, w\Bz\/\B^\ < 30 AU (e.g. 
IAppI et al. 2000). Precession and instabilities therefore lead 
to the fragmentation of the outflow into clumps, which to- 
gether with the generation of episodic ejection, produce a 
series of nested cavities and internal shocks, such as those 
visible in Figure [T^ and f, where faster ejections plough 
through previously launched material. Scaled laboratory as- 
trophysics experiments of compressible, super-magnetosonic 
MHD flows have also shown that although jets may be 
episodic and unstable, multiple ejections of magnetic cav- 
ities and je ts can produce clu mpy flows which remain well 
collimated (|Ciardi et alj|2009h . 



2.3 Outflow suppression and mass accretion. 

We find that the general trend of mass ejection-accretion 
during the collapse of a pre-stellar core, is a decrease in the 
mass ejection efficiency for increasing a, which leads to a 
corresponding increase in the total mass, M core (t), accreted 
onto the central core. Figure [3] shows the total mass in the 
core and outflow (M ou t + M core ) as a function of time. The 
outflow is defined as the plasma with v-R > 0.1 km s~\ 
where R is the unit spherical radius; the core is defined as 
the plasma with a particle number density n > 10 9 cm -3 
and v-R ^ 0.1 km s _1 . From mass conservation the sum 
M ou t + M core is the total infall mass, Mi n j , which at early 



4 Andrea Ciardi and Patrick Hennebelle 




Figure 2. Volume rendering of the outflows at ~ 23000 years, 
for fj, = 5 and a equal to (a) 20°, (b) 45°, (c) 70°, (d) 80° and 
(e) 90° . The colours show increasing speeds (from blue to yellow) 
in an onion-like structure. The black disc-like regions represent 
densities > 10 10 cm' 3 . 

times is seen to be approximately constant, and largely in- 
dependent of both angle a a nd magnetizat ion fi. Its value is 
consistent with that given bv lHunterl(ll977T l. who considered 
similarity solutions of the collapse of an unstable isother- 
mal sphere. The relation between the accreted-ejected mass, 
and the misalignment of rotation axis and magnetic field, 
can be understood by considering the total mass outflow 
rate, f dMout — f pv ■ dS. For steady-state, axisymmetric 
flows it may be written as f dMout ~ f nd& = f kB p -o!S, 
where k — 4npvp /Bp is the mass load per unit time, per 
unit poloidal magnetic flux, and may also be written in 
terms of the density , pa, a t the Alfven surface as k = y/pA 
l|Pelletier fc Pudritzl fl992). We make the simplifying as- 
sumptions that n and BpS ~ B z cosa are constant over 
the surface (S) of the core, so that the integrated mass out- 
flow rate can be approximated by M ou t ~ k$ ~ nB z cos a S. 
We now proceed to establish a scaling relation between the 
mass ejection rate, M ou t, and the mass in the core, M core . 
We assume that the core, of characteristic radius R, is in 
hydrostatic equilibrium, GM core / R oc PR~ 3 , and threaded 
by a constant magnetic flux, <j> oc BR ~ const. The pres- 
sure P and density p scale as P oc p 1 and p oc M core R~ 3 
respectively, and combining these expression one finds: R oc 

M (2- 7 )/(4-3 7)) k M -2/(4-3 7 ) 5 S(xR 2 x M (4-2 7 )/(4-3 7 ) _ 

B oc Afi 2 7- 4)/ < 4 - 37) , and k oc p 1 ' 2 oc M"^ 4 " 37 '. The 
relation obtained for k, is equivalent to assuming that the 
launching speed at the core/disc, v, scales with the Alfven 
speed, B I yf^Kp. From these relations, the mass outflow rate 
is found to scale as a power of the total mass present in 
the core, M ou t oc M^oJe 4 ~ 37 ', which for the adiabatic core 
(7 = 5/3) gives the linear relation: 

Mout = VMcore(t) (1) 

The characteristic accretion-ejection time-scale, n^ 1 = 
T ae = ro/ cos q = const, depends only on the angle a, where 
ro is the value corresponding to the aligned case (a = 0°) 
which is to be determined numerically. We can then com- 
bine Eq. \T\ with the expression for the mass conservation, 
Mi n f = M core + M out , to solve for the mass accreted onto 
the core as a function of time: 

M core (t) = T ae M lnf (l - e - t/T ") (2) 
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Figure 3. Total mass M core + M out for different magnetizations 
p. and angles a. 

and for the total ejected mass: 

M ut(t) = TaeM inf (e~ t/Tae - 1) + M mf t (3) 

Figure [4] shows the simulated infall, core and outflow masses 
for the case p — 5, and for different angles a, together 
with the analytical solutions, which are strictly valid in 
the regime where Mi„f — const. The results show that 
at early times the mass in the core is small, the ejection 
rate is negligible, and mass accretion dominates, M core — 
M inf e- t/Ta ", over the mass ejection, which is given by 
Mout = M in f(l — e -t//Tac ). The model predicts that within a 
time-scale r ae = To/ cos a the core accretes ~ 63% of its "fi- 
nal" mass, which is asymptotically given by Mc = T ae Mi n f 
(for a 7^ 90°). Since r ae increases with the angle a, the initial 
misalignment between the rotation axis and magnetic field 
plays a fundamental role in determining the masses of the 
cores. The limiting case is for a perpendicular magnetic field 
and rotation axis (a — 90°), where no mass is ejected from 
the core, and its growth is at the (time-dependent) mass in- 
fall rate M in j . This case being qualitatively similar to hydro- 
dynamical simulations where no outflows are produced. The 
high ejection efficiency predicted by the model at late times 
is however not energetically favourable. A fraction of the 
outflow fails to become trans- Alfvenic and eventually falls 
back onto the core/disc before being reprocessed. By com- 
bining the mass conservation equation with the definition of 
£, the time when the limiting ejection efficiency occurs (£ = 1 
for jets with negligible enthalpy) can be estimated from the 
relati on: (zui/zu e )^ = 1 — M ou t/ M in f (Fe rreira fc Pelletierl 
1 1995f ) : where w\ ~ 10 — 20 AU corresponds to the core ra- 
dius, and m e ~ 50 — 100 AU corresponds to the outer ra- 
dius of the accretion-ejection structure, where M in f is eval- 
uated. For our model, £ = t [In (■cu e /vji) r ae ] _1 , and the time 
when we expect the the analytical model to overestimate the 
ejected mass is when £ ~ 1, which gives t ~ (1 — 3) x r ae . 
This occurs for the aligned case, which is the only model 
to have reached a high-ejection efficiency, at ~ 24000 years, 
and it is consistent with our estimates. From inspection of 
the velocity fields, it is also evident that part of the out- 
flow is failing to be ejected and falls back as a lower density 
envelope (n < 10 10 cm ~ 3 ). This material is visible in the 
first panel of Figure [4] as an increase in the core's mass with 
respect to the analytical solution. 

The results show that with increasing misalignment, 
the general trend of mass ejection-accretion in collapsing 
pre-stellar cores is to produce higher mass adiabatic cores, 
with lower ejection efficiency (£ oc cos a). Another con- 
sequence is an increase in the bulk speed of the outflow, 
defined as vt — f dVp\v\/ f dVp, with the integral taken 



Outflows and mass accretion in dense cores 5 



over the volume of plasma where v-R > 0.1 km s . Typ- 
ical values at ~ 26000 years are vt(a = 0°) = 1.17 km 
s" 1 , v b (a = 45°) = 1.43 km s" 1 and v b (a = 70°) = 1.76 
km s _1 , with their ratios consistent with the angle depen- 
dence obtained in the analytical model. This increase in the 
bulk speed can be understood by considering the expres- 
sion for the asymptotic je t velocity = Qzoy/ (2 A — 3) 
l|Blandford fe Pavnel (1 19821 ') ). For misaligned, collapsing pre- 
stellar cores two factors contribute to the higher bulk speeds. 
The first is related to the magnetic lever arm A, which de- 
pends on the ejection index through the relation A ~ 1 + ^ 
|Ferreiralll997l '); the second is related to the rotational speed, 
which scales with the core's mass, Qzu ~ y/GM core /R. Thus 
at a given time, the more massive cores obtained for larger 
a, and the corresponding lower ejection efficiency, combine 
to eject faster outflows. Finally we note that following the 
initial launching of the magnetic tower, as £ increases in 
time the outflow from the core/disc is first launched with 
a relatively large magnetic lever arm, and may be thought 
of as magneto-centrifugally driven; at later times, with in- 
creasing £, the lever arm is small and the outflow may be 
better described as magnetic pressure driven. 



3 SUMMARY 

We have presented three-dimensional simulations of the col- 
lapse of pre-stellar dense cores with misaligned magnetic 
field and rotation axis. Imposing non-axisymmetric initial 
conditions has several, important effects on the formation 
of the first adiabatic core and its associated outflows. In 
particular an increasing angle a produces a decrease in the 
ejection efficiency, leading to the suppression of outflows for 
a ~ 90°, and a corresponding increase of the mass of the 
cores. In addition we find that the misalignment of mag- 
netic field and rotation axis produces precessing jets which 
are more susceptible to current-driven instabilities, leading 
to the formation of heterogeneous, clumpy flows. These ef- 
fects are visible even for small a, and indicate that aligned 
configurations may be too idealized, and unable to capture 
the whole dynamics. To interpret the numerical results we 
have developed a model which seems to capture the essen- 
tial features observed in the simulations. In particular we 
find that for low ejection efficiencies the mass ejection rate, 
M ou t, is proportional to the mass of the core, and that the 
reduction of M «t with a is consistent with a decrease of the 
magnetic flux threading the (approximately constant) core's 
surface. 
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